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ABSTRACT 

Despite expanding research activity in gravitational lens modeling, there is no particular software 
which is considered a standard. Much of the gravitational lens modeling software is written by 
individual investigators for their own use. Some gravitational lens modeling software is freely available 
for download but is widely variable with regard to ease of use and quality of documentation. This 
review of 13 software packages was undertaken to provide a single source of information. Gravitational 
lens models are classified as parametric models or non-parametric models, and can be further divided 
into research and educational software. Software used in research includes the GRAVLENS package 
(with both gravlens and lensmodel), Lenstool, LensPerfect, glafic, PixeLens, SimpLens, Lensview, 
and GRALE. In this review, GravLensHD, G-Lens, Gravitational Lensing, lens and MOWGLI are 
categorized as educational programs that are useful for demonstrating various aspects of lensing. Each 
of the 13 software packages is reviewed with regard to software features (installation, documentation, 
files provided, etc.) and lensing features (type of model, input data, output data, etc.) as well as a 
brief review of studies where they have been used. Recent studies have demonstrated the utility of 
strong gravitational lensing data for mass mapping, and suggest increased use of these techniques in 
the future. Coupled with the advent of greatly improved imaging, new approaches to modeling of 
strong gravitational lens systems are needed. This is the first systematic review of strong gravitational 
lens modeling software, providing investigators with a starting point for future software development 
to further advance gravitational lens modeling research. 

Subject headings: strong gravitational lensing, computer simulation, parametric models, non- 
parametric models 



1. INTRODUCTION 

Gravitational lensing has great promise to provide new 
insights into the structure and history of the universe. 
Gravitational lensing has yielded many exciting results 
by mapping dark matter distributions, and the recent use 
of strong gravitational lensing data has added a new di- 
mension to this research (Coe et al. 2010). Gravitational 
lensing is a very active area of investigation, and research 
is highly dependent on computer modeling. In some ar- 
eas of contemporary astrophysics research, there is soft- 
ware that is a de facto standard for many investigators 
(e.g. SExtractor (Bertin and Arnouts 1996) \ GALFIT 
(Peng et al. 2010) 2 , etc.). While gravitational lens mod- 
eling software has been written, there are no standards 
and no easily accessible source of information about ex- 
isting software. The lack of standard software may be 
a virtue of the gravitational lensing community, allow- 
ing more flexibility and creativity. The lack of a single 
standard program makes it more important to compare 
existing software used for modeling strong gravitational 
lenses. Information regarding existing software will be 
helpful to those developing new approaches and inter- 
faces. 

This review was undertaken to identify available strong 
lens modeling software, and review the installation, use, 
and the nature of inputs and data outputs. This paper 
serves as a guide to available software and provides useful 
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1 http : //astroa. physics .metu. edu. tr/MANUALS/ sextractor/ 

2 http : //users . obs . carnegiescience . edu/peng/work/galf it/ 
galf it .html 



information to both new and established investigators in 
this field. The availability of source code may be a use- 
ful starting point for anyone writing their own modeling 
software. 

This paper is organized as follows. In section §2, we 
review the classification of gravitational lens models and 
the methodology used to review the available software. 
In section §3, we review eight software packages that have 
been used extensively in gravitational lens research. Fol- 
lowing this, in section §4, we review five programs that 
are useful in education for gravitational lensing. In sec- 
tion §5 we discuss several factors of importance in se- 
lecting and comparing available software and in section 
§6 we make suggestions for the next generation of soft- 
ware to support future gravitational lens research based 
on this review. 

2. CLASSIFICATION AND REVIEW METHODOLOGY 

Modeling of gravitational lenses starts with a list of 
observables such as relative positions of the components, 
relative fluxes of the images, time delays between the im- 
ages and other lens properties. Gravitational lens model- 
ing can be considered a "forward" or a "reverse" problem. 
In the forward problem, images are predicted based on 
a known source and lensing mass. More commonly, the 
reverse problem is considered, using observed images to 
reconstruct a model of the mass density based on the im- 
ages, usually approached using non-parametric methods. 

In parametric models, a clear physical parameteriza- 
tion is used to construct the model from the outset, while 
non-parametric models are often "grid-based" , although 
there are other methods used. Both parametric and 
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non-parametric models are valuable, since some features 
of individual lenses are model independent while others 
are very model dependent (Saha and Williams 2003). A 
major distinction is whether the calculation is "model- 
based" (parametric) or "model-free" (non-parametric) at 
the start of the process (Coe et al. 2008). A more re- 
cent approach to model classification avoids the confu- 
sion of parametric vs. non-parametric (since all models 
use parameters) and classifies models as light traces mass 
(LTM) or Non-LTM (Coe et al. 2010). 

2.1. Parametric Models 

Parametric simulation models are generally used to 
solve the "forward" problem, taking a source and lensing 
mass, and then predicting the resulting image. Para- 
metric models are also known as simply-parameterized 
models, and assume a physical model which fits the data 
with relatively few defined parameters (Julio et al. 2007). 
In parametric models, the data is fitted to a physical ob- 
ject (e.g. Point Mass, Singular Isothermal Sphere, Singu- 
lar Isothermal ellipsoid, De Vaucouleurs model, Navarro- 
Frank- White Model, etc.) and a model of the lensing 
mass made using that physical object to predict the ef- 
fect on the light from the source. 

2.2. Non-Parametric Models 

Non-parametric models are often used to solve the 
more complex "reverse" problem, also called lens inver- 
sion, of taking a lensed image and from it, predicting the 
nature of the lensing mass. The lens inversion problem 
is complicated by the fact that there are huge degenera- 
cies in the parameter space which make several models 
able to fit observed data (Meneghetti 2012). The de- 
generacy problem is inherent in lens inversion, because 
the constraints on the potential are local (Alard 2009). 
The principle underlying non-parametric gravitational 
lens models is that the effective lens potential and the 
deflection equations are linear functions of the surface 
density (Mcylan ct al. 2006). These models reconstruct 
the mass distribution or lens potential as a map defined 
on a grid of pixels (Julio ct al. 2007). By using a large 
number of parameters these models are very flexible, but 
conversely, the large number of parameters can lead to 
over-fitting the data. The parameters are usually in the 
form of basis functions, and given the large number used, 
the term "non-parametric" is somewhat paradoxical, al- 
though it refers to the lack of a discrete physical model 
being used in the solution of the problem (Liescnborgs 
2009). It is only recently that lens inversion is studied 
using strong gravitational lensing data. 

Lens inversion methods are classified into 3 types: (i) 
model dependent reconstructions, (ii) potential recon- 
struction on a grid and (iii) expansion of the potential 
functions. Degeneracy must be dealt with using any 
of these techniques (Alard 2009). Further examples of 
these techniques include the maximum entropy method 
(Wallington et al. 1996), genetic algorithms (Brewer and 
Lewis 2005; Liesenborgs 2009), Bayesian analysis (Suyu 
et al. 2006; Brewer and Lewis 2006), the semilinear tech- 
nique (Warren and Dye 2003), and perturbative recon- 
struction (Alard 2009). 

The ideal lens inversion algorithm (i) should be free 
of assumptions regarding mass or luminosity distribu- 
tions, (ii) should not depend on prior information, (iii) 



should not produce models that are physically impossi- 
ble, (iv) should be free of uncontrollable parameters, and 
(v) should be extendable to any kind of data (Liesenborgs 
ct al. 2006). The fact that there are so many techniques 
being used to investigate lens inversion, indicates that 
there is no single "best" technique. The theoretical un- 
derpinnings of these various methods have been reviewed 
and shown to be essentially different methods of doing 
the same thing (Brewer and Lewis 2006). 

2.3. Education and Research Modeling Software 

Gravitational lens modeling software can also be ar- 
bitrarily classified into packages used for education and 
those used for research. Those classified for research in 
this review have been used in published studies of gravi- 
tational lensing. Those classified as educational have not 
been used for studies. 

2.4. Review Methodology 

There are two main goals when evaluating gravita- 
tional lens modeling software (Wayth and Webster 2006). 
The first is to determine if the software can recover the 
lens model parameters of a known lens, and with what 
accuracy. The second regards the accuracy of the source 
reconstruction. The accomplishment of these two goals 
must be considered against the background of software 
usability and efficiency. 

Gravitational lens modeling software was identified by 
using search engines on the Internet as well as search- 
ing the literature on arXiv.org. The software reviewed 
here was chosen based on availability for download and 
the ability to install and execute the programs. An- 
other resource for gravitational lens modeling software 
is the Astro-Code Wiki, which includes links for several 
of the modeling codes reviewed here (Astrocode 2013). 
Information regarding release dates and versions was ob- 
tained from the web sites. Each of the software pack- 
ages described here was downloaded, compiled (when 
necessary) and installed. Sample data files were used 
when available, and output shown here is directly from 
the downloaded version. All of the software functioned 
as described by their developers. Documentation was 
also downloaded directly from the internet. This review 
was not intended to be all inclusive; other lens modeling 
packages are not reviewed. Several of these programs, 
although web sites for download were identified on line 
or in publications, were no longer available for download. 
Others had significant issues when attempting to compile 
and execute the programs, precluding their use. There 
are likely other software packages being used by individ- 
uals but not available for download on the internet. The 
software reviewed here is representative of what is avail- 
able in regard to algorithms, types of models available 
and feature sets. Features of the software reviewed re- 
lating to installation and use are summarized in Table 1. 
Features of the software relating to lensing models and 
algorithms are summarized in Table 2. 

3. RESEARCH SOFTWARE 

Strong gravitational lens modeling software for re- 
search use is in this category. The eight software pack- 
ages described here have been used to analyze experi- 
mental data in published studies. 
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forming a calculation requires specifying the details of 
the tiling and the parameters of the mass model. The 
code uses a polar grid centered on the main galaxy. De- 
fault values work with most calculations. Optionally, a 
critical curve grid can be used in which critical curves 
are used to determine where to place radial zones. Vari- 
ous options for plotting the output are available, and the 
output is a file of macros for input into the SuperMongo 
plotting program. The lensmodcl software expands on 
the capabilities of gravlens, to include routines that make 
it easy to fit models to observed lens systems. The cen- 
tral portion of lensmodel computes a x 2 value for a set of 
models. A wide variety of lens models are available in the 
catalog (Kccton 2001a). The software can also calculate 
the Hubble constant using time delays. The available 
tutorials are very helpful for illustrating and explaining 
the complex commands that are used in gravlens and 
lensmodel. 

Gravlens was used to model strong gravitational lenses 
in an interesting study which examined the cusp and fold 
relations as a gauge of substructure ( Aazami and Natara- 
jan 2006). This was a theoretical investigation to look 
at the sensitivity of these relations to the presence of 
substructure in the lens. The authors found that the 
fold relation is a more robust indicator of substructure 
than the cusp. Gravlens provided an excellent analyt- 
ical tool for this study. A newly discovered elliptical 
galaxy at z=0.0345 was investigated and modeled with 
gravlens/lensmodel (Smith et al. 2005). The model was 
created with a SIE, and was found to be in close agree- 
ment with the light distribution observed. 

In a study of a bright strongly lensed z=2 galaxy in 
the SDSS DR5, Lin and colleagues used Lensvicw to 
model a lensing system (Lin et al. 2009). However, since 
Lensview uses the full image information, they employed 
gravlens/lensmodel to fit an SIE model using only the im- 
age positions. The resulting model showed a very good 
fit to the image positions. This study is one of the very 
few direct comparisons of strong gravitational lens mod- 
eling software in the literature. 

The ten image radio lens, B 1933+503 was modeled us- 
ing gravlens by Cohn and coworkers (Cohn et al. 2001). 
The mass distribution of this system was modeled with 
a wide variety of parameterized ellipsoidal density dis- 
tributions. The models were constrained using the rel- 
ative positions of the lens galaxy and the lensed images 
and the flux ratios between the images. The mass dis- 
tribution was concluded to have an approximately flat 
rotation curve based on this analysis. The gravlens pro- 
gram was used to model the cluster lens MS0451. 6-0305 
by Alba and colleagues (Berciano Alba et al. 2007). A 
simple elliptical lens model (SIE) with external shear was 
used. They used an NFW profile for the cluster mass dis- 
tribution which was consistent with observations. The 
model reproduced the positions of the Extremely Red 
Object images and the radio images very well. This 
model was used to test if the configuration of the ob- 
served radio emission could be understood as a result of 
gravitational lensing, and was able to explain the mor- 
phology of the radio map as a result of the three lensed 
background sources. Gravlens was also used to model a 
strongly lensed Lyman Break galaxy at z=2.73, identi- 
fied in the SDSS DR4 imaging data (Allam et al. 2007). 
The authors assumed a SIE and used gravlens/lensmodel 



GravLensHD 2011 No Yes iOS 

G-Lcns 1998 No Yes DOS 

Gravitational Lensing 2002 No Yes Win/HP-49 

lens 2002 Yes No MATLAB 

MOWGLI 2013 No Yes Java 

TABLE 1 

A SUMMARY OF SOFTWARE FEATURES OF GRAVITATIONAL LENS 
MODELING SOFTWARE. SOURCE INDICATES THAT THE SOURCE CODE 
IS PROVIDED, Exec INDICATES THAT A DOWNLOADABLE 

executable file is provided and platform indicates the 
computing platform that the executable works on. 
(PPC=PowerPC, Dates are the current version available, 
exec— available executable code) 

3.1. gravlens / lensmodel 

3.1.1. Software features 

The GRAVLENS- Software for gravitational lensing 
package was first released in 2001, and is now in version 
1.99 dated 2008 3 . The software is available as a down- 
load, with executable files provided for PPC and Linux 
architectures. Source code is not available. Documenta- 
tion is also available as a 101 page user manual (Keeton 
2004). The lensmodel program is available as a down- 
load, with executable files provided for PPC and Linux 
architectures. Source code is not available. The lens- 
model software uses the gravlens kernel and adds func- 
tionality, and was also described by Keeton (2001b). A 
tutorial is provided with the user manual. The user in- 
terface is character based and input is through a text 
file with text commands. There are two on-line tutori- 
als that illustrate many of the features of the software. 
While gravlens includes basic lensing calculations, lens- 
model includes added routines to model strong lenses 
(Keeton 2004). 

3.1.2. Lensing features 

The lensmodel software uses a parametric model, and 
has been used in a number of research studies (Cohn 
et al. 2001; Allam et al. 2007; Berciano Alba et al. 2007; 
Lin et al. 2009; Smith et al. 2005; Aazami and Natara- 
jan 2006). GRAVLENS was first described along with a 
new algorithm that allowed calculation of lensing param- 
eters for a generalized mass model, allowing implemen- 
tation of multiple parameterized models (Keeton 2001b). 
GRAVLENS is a sophisticated software package that is 
accompanied by a large catalog of parametric lensing 
models (Keeton 2001a). Both programs in GRAVLENS 
includes a wide range of basic lensing calculations which 
are based on tiling of the image. The heart of the code 
is a general algorithm for solving the lens equation and 
is fully described in Keeton (2001b). This algorithm in- 
volves tiling the image and source planes, and using these 
tiles to determine the number and approximate positions 
of all lensed images associated with a given source. Per- 

3 http : //redf ive . rutgers . edu/~keeton/gravlens/ 
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to perform fits to the data. The fitted values showed ex- 
cellent agreement with observed values from the SDSS 
DR4 database. 

3.1.3. Summary 

Both gravlens and lensmodel require no installation as 
they are provided as executable binary files. They are 
somewhat complicated to use because of the character 
based interface and extensive command set. They would 
benefit from a graphical user interface which is described 
in a web site but the code has not been available (Al- 
faro 2008) 4 . The lack of source code is also a limit- 
ing factor for those wishing to study the computational 
techniques used. Finally, the programs are limited by 
requiring an outside plotting package to view graphical 
results. The required plotting package is expensive and 
somewhat outdated. Both gravlens and lensmodel have 
an extensive catalog of available mass models that can 
work with complex datasets, and perform very well in 
regard to comparing the models generated with observa- 
tional results. Overall, they are excellent strong gravita- 
tional lens modeling software packages. 

3.2. Lenstool 

3.2.1. Software features 

Lenstool was first written in 1993, and is now released 
as version 6.7.1, dated 2006 5 . Lenstool is distributed 
as source code with dependencies on WCSTOOLS, PG- 
PLOT and CFITSIO, and the latest version is available 
on the Lenstool Project web site (Kneib 2012). Installa- 
tion is somewhat tedious but can be accomplished with 
available standard libraries, and can be installed in Linux 
or OS X. The user interface is character based and in- 
put is through a text file with text commands. A 61- 
page user manual as well as a 41 page document en- 
titled "Lenstool for Dummies" are available for down- 
load (McCourt 2006). The documentation is excellent 
for this complex and comprehensive strong gravitational 
lens modeling software. 

3.2.2. Lensmg features 

Input to Lenstool is a character based input file. Each 
line consists of a command and appropriate data. Key- 
words are either first identifiers or second identifiers. A 
number of input files are necessary. The first is a PAR 
file, containing the basic parameters for the model, a list 
of the arcs for which Lenstool will predict counter im- 
ages, and specification of requested optimization. The 
remainder of the input includes multiple image files and 
a cluster members file. The available options and com- 
mands are extensive, and testing the software with the 
files available for download is an excellent way to gain 
familiarity with the complex input required. 

Lenstool is referenced in at least 11 manuscripts (e.g. 
Julio ct al. 2007; Kneib et al. 1996; Richard et al. 2008; 
Limousin et al. 2007; Richard et al. 2009; Julio and Kneib 
2009). These manuscripts are listed on the Lenstool 
website (Kneib 2012), and many of them include down- 
loadable files of the lens models used in the research. 

4 http : //cinespa.ucr . ac . cr/sof tware/xf gl/index .html 

5 http : //www . oamp . f r/cosmology/lenstool/ 



Strong lens galaxy clusters are modeled with paramet- 
ric methods, and ranked using Bayesian evidence. Al- 
though Lenstool was initially developed in 1993 with 
a downhill x 2 minimization, modeling of complex sys- 
tems become inefficient due to the sensitivity of the tech- 
nique to local minima. The computational method used 
by Lenstool was then changed to use a publicly avail- 
able Markov Chain Monte-Carlo sampler, avoiding local 
minima in the likelihood functions. The merits of this 
method on simulated strong lensing clusters is demon- 
strated by Julio et al Julio ct al. (2007). 

Using a multi-scale model with a hybrid approach of 
LTM and non-LTM modeling in the Lenstool software, 
Julio and Kneib were able to model Abell 1689, but only 
for a limited subset of images (Julio and Kneib 2009). 
The key feature is that Lenstool uses a multi-scale model, 
allowing sharper contrast in regions of higher density. 
This arrangement of potentials of different sizes allows 
Lenstool to produce a high-resolution model with a min- 
imum number of parameters. In this study, a mass recon- 
struction was created using Lenstool. The model com- 
bined a grid of radial basis functions and galaxy scale 
clumps with cluster member galaxies. A grid was built 
from a mass map based on 2 cluster-scale and 60 galaxy- 
scale clumps of mass, instead of the 190 galaxy-scale 
clumps used in previous studies (Limousin et al. 2007). 
A catalog of 28 images in 12 systems of multiple images 
were selected for this analysis. There were 122 parame- 
ters, which took 15 days to produce 2000 MCMC samples 
on a 2.4GHz processor. The results of this study con- 
firmed the ability of a multi-scale model to be used as 
a lens model. The authors report that the errors be- 
tween the positions of observed and predicted images 
were halved, compared to previous studies. 

The nature of the mass distribution in Abell 1703 
was studied using Lenstool (Richard ct al. 2009). This 
demonstrated the ability to model the inner mass dis- 
tribution of massive galaxy clusters. This study used 
a spectroscopic survey to confirm photometric redshifts 
and precisely constrain the mass distribution in Abell 
1703. Lenstool was used to constrain a parametric mass 
model with the identified multiple systems. The posi- 
tions of the multiply imaged systems were used to op- 
timize parameters describing the mass distribution, us- 
ing a Pseudo Isothermal Elliptical Mass Distribution 
(PIEMD) using the profile derived from photometry. The 
Bayesian approach in Lenstool provides a large number 
of models which sample the probability density function 
of all the parameters. Results were compared to previ- 
ous weak- lensing studies, and were found to have a very 
good fit. This strong lensing analysis using Lenstool and 
a simple NFW component for the large scale dark matter 
distribution, was able to reproduce the large number of 
images in Abell 1703, as well as demonstrate consistency 
with previous weak-lensing analyses. 

3.2.3. Summary 

Lenstool is a comprehensive program for gravitational 
lens modeling. The software is relatively easy to install 
and use, and is accompanied by extensive documenta- 
tion. It has been extensively used in the modeling of 
observed gravitational lenses. The availability of the in- 
put files from previous studies increases the utility of this 
software. The software uses a novel technique to combine 
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the strengths of both LTM and non-LTM models, by us- 
ing a multi-scale model that allows sharper contrast in 
areas of high density. The flexibility of this approach al- 
lows improved prediction of image positions compared to 
previous studies. 

3.3. LensPerfect 

3.3.1. Software features 

LensPerfect was released in 2007 and is available for 
download 6 . The software is written in Python and the 
source code is included. Dependencies include Numpy, 
Scipy, Matplotlib and Pyfits. Installation of the software 
is slightly complicated by the fact that the software will 
not work with the latest version of Python. There are 
numerous sample data files provided to demonstrate the 
features of the software. Output from analysis of one 
of the sample files is shown in Fig 1. Documentation is 
available on two web sites including one which details 
installation and use of the software (Coe 2012), although 
there is no separate user manual. The user interface is 
character based and input is through a text file with text 
commands. 

3.3.2. Lensmg features 

The software was developed by Coe (Coe et al. 2008) 
and has also been used in one more research study (Coe 
ct al. 2010). LensPerfect uses a parametric model but 
is also "model-free" as described by its developers, who 
further characterize it as non-LTM. LensPerfect solutions 
are given as sums of basis functions. While most para- 
metric models are based on a physical object, the basis 
functions used by LensPerfect have no physical interpre- 
tation. Input to the program is via a text file and graph- 
ical output is shown immediately on the display. 

LensPerfect represents a new approach to gravitational 
lens mass map reconstruction, and is the first method 
to do so using strong gravitational lensing data (multi- 
ple images). This new approach uses direct inversion to 
obtain assumption-free mass map solutions which per- 
fectly reproduce all multiple image positions. This was 
developed using a new mathematical approach, using a 
curl-free interpolation of vectors given at scattered data 
points. One of the key features in any model is the mea- 
sure of " physicality" of the model. The developers of 
LensPerfect use a new measure of physicality, with the 
following traits: (i) positive mass everywhere within the 
convex hull, (ii) low mass scatter in each radial bin, (iii) 
no "tunnels", (iv) overall smoothness and (v) average 
mass in radial bins decreases outwards. The only rigid 
constraint among these is the first trait, requiring posi- 
tive mass. A complete discussion of these traits is avail- 
able in the reference (Coe et al. 2008). 

LensPerfect provides an accurate mass map even when 
there are many lensed galaxies, by using several novel 
approaches. A weighted average of predicted source po- 
sitions is used to determine each new source position. 
The solution is rebuilt at each iteration as new sources 
are added. This process is fast, and results in an accurate 
mass map. Both source and image positions are always 
perfectly constrained. 

6 http : //www. its . caltech . edu/~coe/LensPerf ect/download/ 
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Fig. 1. — Sample output from LensPerfect using a supplied test 
dataset 

The galaxy cluster Abell 1689 is one of the most stud- 
ied gravitational lens systems, and thus is ideal for com- 
parisons among lensing models. The large number of 
multiple images in Abell 1689 also make detailed anal- 
ysis a challenging problem. In a followup investigation 
after introducing the method and software, LensPerfect 
was applied to the analysis of Abell 1689 using the po- 
sitions of 168 multiple images. The non-LTM models 
from LensPerfect were able to reproduce the observed 
input positions of 168 multiple images of 55 knots resid- 
ing within 135 images of 42 galaxies (Coe et al. 2010). 
The computing problem associated with analyzing Abell 
1689 is obvious, since it has 100+ strong lensing features. 
The software must produce a mass model with correct 
amounts of mass to deflect light from 30+ background 
galaxies into multiple paths to arrive at the 100+ ob- 
served positions. LensPerfect did this, using direct ma- 
trix inversion to find solutions based on the input data. 
The optimization took two weeks to run on a Macbook 
Pro (Apple Corp., Cupertino CA) laptop computer. This 
study demonstrates the robustness of the algorithm and 
its computing efficiency. 

A mass map of Abell 1689 produced by LensPerfect 
using NFW and Sersic fits had a recovered mass pro- 
file which matched the input mass profile extremely well 
(Coe et al. 2010). The NFW fit parameters compare very 
well to previous studies of Abell 1689 using strong lens- 
ing data as well as studies using a combination of strong 
and weak lensing data. 

3.3.3. Summary 

LensPerfect is a full-featured program that is extremely 
capable of producing a mass map based on lensed images, 
even when there are many galaxies as demonstrated in 
the followup study of Abell 1689. It is distributed as 
source code, written in Python, which is available for 
many platforms. The new approach taken by LensPer- 
fect results in fast and accurate modeling, even of com- 
plex systems. The output is graphical and displayed 
directly. LensPerfect does an excellent job of accurate 
reconstruction of the sources and lensed images. The 
ability of LensPerfect to accurately produce a mass map 
using strong lensing data represents an important ad- 
vance in gravitational lens research, and will likely lead 
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to further advances. 

3.4. glafic 

3.4.1. Software features 

The current release of glafic was written in 2012, and 
is provided as a downloadable executable file which is 
unpacked from a .tar file 7 . Several sample input files 
(also in a single .tar file) and a detailed 51 page Users 
manual are provided. Source code is not provided. There 
is no installation procedure since it is provided as an ex- 
ecutable file for OS X. The program runs without any 
modifications necessary. The user interface is character 
based and input is through a text file with text com- 
mands. The program runs in the command line interface 
using OS X. There is no graphical interface. 

3.4.2. Lensing features 

Glafic uses a parametric model that can be used for 
a wide variety of gravitational lensing analyses (Oguri 
2010). It includes computation of lensed images for both 
point and extended sources, handling of multiple sources, 
a wide variety of lens potentials and a technique for mass 
modeling. Commands are entered in a simple text file, 
which begins with a list of primary parameters (omega, 
lambda, Hubble, lens redshift z, pixel size, etc.) and then 
an optional list of secondary parameters (optional data 
files, output format desired, extended source normaliza- 
tions, etc). Point sources are defined simply by their 
redshift and x- and y-coordinates. Each lens is defined 
by the lens model and seven parameters. A large catalog 
of lens models is available (including point mass, Hern- 
quist, NFW, Einsato, Sersic, etc.). Extended sources can 
be Gaussian, Sersic, top hat, Moffat or multiple sources. 
After defining the parameters and the lens models, pa- 
rameters to be varied in the x 2 minimizations are spec- 
ified. Following this, the desired commands are issued 
such as computing various lensing properties, Einstein 
radius, write lensing properties to a FITS file (see sample 
output in Fig 2), etc. Various types of optimization are 
permitted. The commands are well described and illus- 
trated in the User's manual. Commands can be entered 
as a batch using an input file, or entered on a command 
line. A number of sample data files are provided which 
illustrate a number of the major features of glafic. 

Glafic was used to perform a strong lens analysis of 
SDSS J1004+4112 (Oguri 2010). This is a particularly 
interesting quasar lens because it is one of only two 
known examples of a cluster-scale quasar lens, and con- 
tains multiply imaged galaxies at z ~ 3. The authors 
include an indirect comparison with multiple previous 
mass models of this interesting lens. This study used 
a parametric model, with the main halo of the lensing 
cluster modeled with the generalized NFW profile. A 
standard \ 2 minimization was used to find the best-fit 
mass model. The best-fit radial mass profile generated 
is in good agreement with strong lensing data inferred 
from Chandra X-ray observations. The model used sev- 
eral new constraints including positions of spectroscop- 
ically confirmed multiple imaged galaxies, time delays 
between quasar images, and faint central images. The 
model thus generated was able to successfully reproduce 
all observations including time delays. 

7 http : //www . slac . stanford.edu/-oguri/glafic/ 
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Fig. 2. — FITS image output by glafic for a sample model with 
a single point source, two lenses and two extended sources 
3.4.3. Summary 

The glafic program is a well-designed, flexible and easy 
to use lens modeling program. It includes a wide vari- 
ety of lens models and is extremely flexible in the types 
of models it accepts. The User's manual is helpful and 
complete. Output consists of a variety of possible data 
files including FITS images. Glafic can accurately re- 
cover lens model parameters of known lensing systems. 

3.5. PixeLens 

3.5.1. Software features 

PixeLens was written in 2007, and is available to run 
in a web browser or as a stand alone java applet which 
can be downloaded from the website 8 . There is no in- 
stallation needed. This program was used at the AN- 
GLES School for Gravitational Lensing in 2007. The 
source code is provided in the downloadable .jar file. In- 
formation is also available on the website. Sample input 
data is provided. A paper (19 pages) and tutorial (14 
pages) detailing the software and underlying theory are 
available on the website, and serve as excellent documen- 
tation with numerous examples. The user interface is a 
single window with several entry panels in the window. 
Input is done in the window, or through an input file. 

3.5.2. Lensing Features 

This program uses a non-parametric model, with 
Bayesian statistics. Input to the program consists of 
model constants (red shifts, pixel size, etc) and image 
data. The radius of the mass map in pixels and the red- 
shifts for the lens and source must be given. Optional 
inputs include setting the mass map as symmetric or 
asymmetric, the radius of the mass map, external shear, 
the number of models allowed, and several others. Im- 
age data includes x- and y- locations and the time delays. 
Type of output can be selected as text or eps files. Sam- 
ple output is shown in Fig. 3, showing the image directly 
on the display. 

Saha et al. (2006b) used PixeLens using non-LTM 
models to perfectly reproduce some of the data for 
Abelll689, but computational limitations restricted Pix- 
eLens to fitting 30 multiple images at a time, because of 

8 http : //www . qgd.uzh. ch/projects/pixelens/ 
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Fig. 3. — Sample output from PixeLens 

the requirement for exact fits to the image data. This 
could be considered a virtue allowing the data to be split 
into two sets as pointed out in the study. In this study, 
PixeLens had been enhanced to use multiple-source red- 
shifts. PixeLens had some advantages in the models gen- 
erated, particularly in handling the inherent problem of 
degeneracy. While some software uses one or a few mod- 
els, PixeLens generates a large ensemble of models which 
explore the possible mass distributions that can repro- 
duce the data. Data from PixeLens was compared to in- 
dependent datasets for consistency (Sana et al. 2006b). 
PixeLens has also been used to model the giant quadru- 
ple quasar SDSS J1004+4112 (Williams and Saha 2004; 
Saha et al. 2006b). In addition to the lensing data (image 
positions) , six kinds of constraints were applied to limit 
the ensemble to lenses that could plausibly be galaxies 
or clusters. The result was free-form reconstructions al- 
lowed detection of structure in the lens associated with 
cluster galaxies (Williams and Saha 2004). J1004+4112 
was reconstructed using 13 images from 4 sources. 

3.5.3. Summary 

PixeLens is an easy to use program for non-parametric 
modeling of gravitational lenses. There is ample doc- 
umentation available. It is based on reconstruction of 
a pixellated mass map by generating large ensembles of 
models with a Metropolis algorithm. It can model several 
lenses simultaneously, which is a rare feature of modeling 
software. The code has been tested with a number of fake 
models and correctly recovers both the mass distribu- 
tions of the lens and HO within uncertainties. PixeLens 
has been used in several published studies to analyze 
gravitational lenses (Saha et al. 2006a, b). 

3.6. Simp Lens 
3.6.1. Software features 

SimpLens was written in 2003, and is available to run 
in a web browser or as a stand alone java applet which 
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Fig. 4. — Sample data entry and output from SimpLens 

can be downloaded from the website 9 . There is no in- 
stallation procedure needed. The source code is provided 
in the downloadable .jar file. The user interface is a sin- 
gle window which displays graphs of mass, potential and 
arrival time. There is no formal written documentation, 
but the web site provides adequate information. 

3.6.2. Lensing features 

SimpLens uses a non-parametric model, and the algo- 
rithm is explained in an article written by the developer 
(Saha and Williams 2003). That article also demon- 
strates some of the analytic capabilities of SimpLens. 
The user can enter values for a,b,h,gl, g2, n, eps, r, x, 
and y. The parameters a, b, and h specify an elliptical 
mass distribution, while gl and g2 refer to the external 
mass. By clicking the computer mouse, the source posi- 
tion can be changed (see Fig 4 left panel) and the effect 
of that change instantly visualized on the display. The 
graphs display caustics, saddle point contours and criti- 
cal curves. A number of sample data sets are provided 
which are easily selected in a pull-down menu. The vari- 
ous parameters from the sample data can then be varied 
to instantly observe the effect of the alteration on the 
three graphs. For example, the sample data allows one 
to easily observe the effect of added shear. 

3.6.3. Summary 

SimpLens is a simple interactive program which allows 
one to instantly demonstrate the effect of changing pa- 
rameters on caustics, saddle points and critical curves for 
mass, potential and arrival time. It is very easy to use 
and very instructive. 

3.7. Lensview 
3.7.1. Software features 

Lensview was first released in 2006 and he current ver- 
sion 1.1.2 is dated 2008 10 . This software is distributed 
as source code, and has dependencies on CFITSIO, GSL 
and the FFTW libraries, as well as an optional depen- 
dency on fastscll. It can be installed in Linux or OS 
X. There is a web site where the software can be down- 
loaded, with information regarding installation and use, 
but there is no separate user manual (Wayth 2012). The 
user interface is character based and input is through a 
text file with text commands. The software is described 

9 http : //www .physik. uzh. ch/~ps aha/lens/ simplens . php 

10 http : //cira. ivec . org/dokuwiki/doku.php/staf f /rwayth/ 
lensview 
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as running in two modes, "simple projection" and "nor- 
mal" (fit an observed image). 

3.7.2. Lensmg features 

Lensview is based on the LensMEM algorithm 
(Wallington et al. 1996), which finds the best fitting 
lens mass model and source brightness distribution us- 
ing a maximum entropy constraint (Wayth and Webster 
2006). Lensview is used to study lens inversion, obtain- 
ing a model of the source based on lensed images. It 
uses a parametric model for solving the lens inversion 
problem. The main features of the software include: (i) 
projection between image and source plane conserving 
surface brightness, (ii) compound lens models which can 
contain several basic components, (iii) reconstruction of 
the unlensed source brightness profile and corresponding 
model image using a non-parametric source, (iv) itera- 
tive source reconstruction process incorporating a maxi- 
mum entropy metric, and (v) statistical evaluation of the 
model image given by the current lens model parameters. 

In simple mode, Lensview takes as input: a source, 
a lens model, a data image (optional) and projects the 
source onto the image. The projected image and x 2 
are output. A log file is also produced. A variety of 
lens model components are supported, including pseudo- 
isothermal potential (PIEP), SIE, point mass, power-law 
mass distribution, NFW, constant density mass sheet, 
external shear, exponential disc and a Sersic mass distri- 
bution. Sample templates for each of these mass distri- 
butions are given on the website (Wayth 2012). To use 
the software, one needs a data image, a PSF image and a 
mask image, all of which need to be the same size. There 
are no sample data files provided with the software, but 
there are some examples available through the website. 
Output from the program can be a FITS image, or a file 
for input into IDL and plotting the \ 2 surfaces. 

In their study, Wayth and Webster first used Lensview 
with simulated optical images to determine optimal pixel 
size, accuracy of lens model parameters, distinguishing 
power of the image and overall goodness-of-fit (Wayth 
and Webster 2006). Model images were created with 
a two-component source having two regions of different 
size and peak brightness. A PSF was generated and con- 
volved with the simulated image. They found that an 
image to source pixel scale ratio of at least 2 is required 
to reproduce the data and a ratio of 3 was optimal. A suf- 
ficiently large source plane of approximately 15x15 pixels 
was adequate. All lens model parameters were well re- 
covered using this data. A comparison of various models 
showed that the image could distinguish between lens 
models although the differences were quite small. 

Lensview was used to analyze the Einstein Ring ER 
0047-2808 (Wayth et al. 2005). This study modeled 
the system using six different models, including PIEP, 
SIS+7, SIE, SPEMD, M/L and NFW. This study showed 
that the SIS+7, M/L and NFW (25kpc) could not repro- 
duce the data, but that the other models could. The dif- 
ferences between the other, successful models, were sub- 
tle. The data showed that the lensed image has four dis- 
tinct bright regions, and that Lensview was able to gen- 
erate a best-fitting image and reconstructed the source 
brightness profile using a non-parametric model. 

Lensview was also used to analyze the optically lensed 
arc HST J15433+5352 (Wayth and Webster 2006). This 



lens was modeled with a PIEP model, including an ex- 
ternal shear. The final source plane used was 10x10 pix- 
els, with a source-to-image plane pixel scale ratio of 1/2. 
The resulting model was somewhat surprising in that a 
purely elliptical model reproduced the data equally to 
a model including shear. They found a critical radius 
(b) of 0.525, which differed somewhat from data previ- 
ously reported with a critical radius of 0.58 (Knudson 
et al. 2001). Some of this difference could be accounted 
for the fact that the Lensview study is a non-parametric 
model, while Knudson et al Knudson ct al. (2001) used 
a parametric model. 

Lin et al used Lensview to study a bright strongly 
lensed galaxy in the SDSS DR5 (Lin et al. 2009). 
Lensview provided an excellent model, but as the au- 
thors pointed out, Lensview uses the full image informa- 
tion which precludes determination of how well the image 
positions are determined. They then used Lensmodcl to 
fit an SIE model using only the image positions. This 
study illustrates both the strengths and weaknesses of 
this approach, but more importantly also suggests the 
need for using multiple approaches to modeling. 

3.7.3. Summary 

Lensview is a comprehensive modeling program, with 
a large number of features and options. It is not straight- 
forward to use, although the sample data files and sug- 
gestions provided on the website do facilitate gaining pro- 
ficiency with the software. Due to its comprehensive na- 
ture, it is possible to specify very complicated lens models 
based on one or more components. 

3.8. GRALE 

3.8.1. Software features 

The current version of GRALE is 0.9.0, released in 
2008 11 . It is provided as source code which must be com- 
piled and linked using the CMAKE utility. The installa- 
tion procedure is fairly straightforward, and the GSL and 
CFITSIO libraries are necessary. The software can be in- 
stalled in Linux, or OS X. The GRALE library is easily 
run under GRALESHELL, which provides an interactive 
environment. The software runs in a single window and 
text commands are entered in a panel at the bottom of 
the window. Output is obtained using GNUPLOT (see 
Fig. 5) or optionally as a FITS image. Documentation 
consists of a website dedicated to this software, but there 
is no separate user manual. 

3.8.2. Lensing features 

GRALE can be used to simulate gravitational lenses 
and to invert lensing systems. The GRALE algorithm 
uses a non-parametric technique to infer the mass distri- 
bution of a gravitational lens system with multiple strong 
lensed systems (Licsenborgs ct al. 2006). To start sim- 
ulating a gravitational lens, the user first decides which 
type of lens to use. There are a variety to select from in- 
cluding a point mass lens, a SIS, SIE, projected Plummer 
sphere, square shaped region of constant density, a two 
dimensional gaussian density profile, and others. The 
lens parameters include mass, distance and others as ap- 
propriate. Commands and parameters are entered line 

11 http: //research. edm.uhasselt .be/~j or i/page/index .php? 
n=Physics . Grale 
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by line in GRALESHELL. Next, the distances to the 
source plane are specified and then the mapping from 
the image plane to the source plane is specified. Sources 
can then be added. A sample lens model is shown in 
Fig. 5. Lens parameters can be optionally saved and 
retrieved. The program directly outputs a file which is 
used by GNUPLOT to generate visual output. 

GRALE has been used in a number of published stud- 
ies to analyze various lenses. GRALE functions well in 
modeling systems with few images, and thus less infor- 
mation (Liesenborgs et al. 2007). GRALE was used in 
a study using strong lens modeling to search for dark 
matter (Liesenborgs et al. 2008b). 

GRALE was used to model the system CI 0024+1654, 
and infer the mass maps using a non-parametric tech- 
nique (Liesenborgs et al. 2008a). This represents the first 
time that image information alone was used to recon- 
struct the mass distribution of this cluster in the strong 
lensing region. No information about the positions of 
cluster members was used. Image data for sources A and 
B as previously described was used for the inversion pro- 
cedure to reconstruct the source. Source A was mapped 
onto five images and source B mapped onto two images. 
A mass map was obtained by averaging 28 solutions. The 
only bias in the technique is that the user must specify a 
square shaped area for the algorithm to search mass dis- 
tribution, assuming that no mass is outside the bound- 
aries of that region. Using the inversion procedure in 
GRALE, an averaged mass map was obtained that dis- 
played the features seen in the ACS images. 

Models of the well known system SDSS J1044+4112 
were also created using GRALE (Liesenborgs et al. 2009) . 
This study looked at five spectroscopically confirmed im- 
ages and position information was used from existing 
studies. One of the images was uncertain (A5) and was 
not included in the first inversion. A second inversion 
was then performed with the addition of image A5, which 
allowed reconstruction of the source shapes after project- 
ing back onto the source planes. Calculation of the total 
mass within 60Kpc and llOKpc compared well to results 
in existing studies. Comparison was also made to the 
best fit NFW model, and the configuration of the cluster 
corresponded to that reported previously (Oguri et al. 
2004). 

3.8.3. Summary 

GRALE is easily installed and straightforward to use. 
Commands are text-based and relatively intuitive. A va- 
riety of lens models is available for simulations. A large 
number of lens and source parameters can be entered 
by the user, allowing one to perform complex lens simu- 
lations. In several studies, GRALE did an excellent job 
of recovering lens parameters and accurately reconstruct- 
ing the source parameters. GRALE has been tested with 
models of C 0024+1654 and SDSS J1004+4112 and gives 
results compatible with previous studies. 

4. EDUCATIONAL SOFTWARE 

Programs suited for education are useful to demon- 
strate the basic principles of lensing. Categorization as 
educational software is arbitrary and based on the fact 
that there were no studies found in the literature that 
used these four software packages for data analysis in 
published studies. 
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Fig. 5. — Output from GRALE for one of the supplied sample 
data files 

4.1. GravLens HD 

4.1.1. Software features 

GravLensHD is available for iOS devices and is ob- 
tained from the Apple (Cupertino CA) App Store, at no 
cost 12 . The source code is not provided. There is a help- 
screen with some suggested exercises to demonstrate the 
features of a gravitational lens. The web site has some 
additional information (Rykoff 2012). There is no for- 
mal documentation. The interface is limited to a single 
touch screen. Installation is simple using the App Store 
application. 

4.1.2. Lensing Features 

This software uses a single mass as a lens, and is cat- 
egorized as a parametric model. The software uses a 
simulated source and a foreground lensing mass to draw 
arcs or an Einstein ring as one moves the lensing mass 
relative to the star digitally on the touch screen. The 
lensing mass can be made invisible leaving the lensed 
image of the source. The ellipticity of the lensing mass 
can be changed through a few preset shapes. The size 
of the lensing mass can be changed with on-screen ges- 
tures. The images can be saved to the devices library.The 
background image can be changed to any picture of the 
user's choice, giving an excellent way to show the effect 
of a lens on a known image. Some technical material is 
also available within the application. It uses a singular 
isothermal sphere (SIS) model with external shear. The 
maximum shear permitted is 0.6. 

4.1.3. Summary 

In summary, this software provides an excellent demon- 
stration of the effect of a gravitational lens, and may 
serve as a tool to stimulate young minds to think about 
astronomy. The concepts used are the same as those for 
any gravitational lensing system. There is no option for 

12 Available on iOS devices though the App Store (Apple Co. 
Cupertino CA) application 
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numerical input, and all interaction with the software is 
through the touch screen. Further information is avail- 
able on a web site dedicated to this software (Rykoff 
2012). 

4.2. G-Lens 

4.2.1. Software features 

G-Lens was written by Boughen in 1998 (as part of an 
undergraduate honors thesis), and is available as an on- 
line download as a single executable file 13 . The program 
will only run in a DOS or early Windows environment, 
but can be used with more modern operating systems in 
DOS emulators, such as DOSBox which runs in OS X 
(Veenstra 2012). Source code is not provided. There is 
no installation, as one needs only to execute the .exe file 
provided. The software was reviewed in the June 2000 
issue of Sky and Telescope (Boughen 2000) . Along with 
the software, there is a one-page instruction sheet avail- 
able, but no other documentation. The user interface 
is a single character based text box where the values of 
parameters are entered. 

4.2.2. Lensing Features 

The lens mass is varied and is modeled as a single point 
mass, characterizing this software as a parametric model. 
Input to the program is performed by providing the lens- 
ing mass (in solar masses), the lens and source distances 
(in Mpc) and then selecting among three options for ge- 
ometry of the lens (circle, ellipse and grid). Output is 
immediately shown on the display. The image can be 
optionally printed. 

4.2.3. Summary 

This software simulates a simple point mass lens and 
allows numerical input regarding basic system geometry. 
The output is graphical showing the image resulting from 
the lens designed by the user. The ability to specify basic 
parameters numerically makes this a good program for 
demonstration of basic lens effects. 

4.3. Gravitational Lensing 

4.3.1. Software features 

Gravitational Lensing was written in 2002 and is dis- 
tributed as an executable file which runs in Windows, 
and also as a C program for use on an HP-49 calculator. 
The software is available for download on the website 14 . 
There is no installation, as one needs only to run the 
.exe file provided. Source code for the Windows version 
is not distributed. Documentation is available to explain 
the science of lensing, written in Dutch. There is no 
documentation for running the software. The interface 
is a single window, and a "File" menu. The program 
has an "About" screen but no on-line help. The options 
in the "File" menu allows changing parameters of the 
lens model. After changing the parameters, the result is 
displayed in the window. 

4.3.2. Lensing features 

13 http : //astrotips . com/sof tware/g- lens 

14 http : //www. kwakkelf lap . com/gravlens . html 
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Fig. 6. — Sample output from Gravitational Lensing 

The mass can be varied on the data entry screen and is 
modeled as a point mass, characterizing this as a para- 
metric model. Sample output is shown in Fig 6. The 
mass of the source and lens can be set individually (as 
multiples of solar masses). Distances from the observer 
to source and lens are set in Mpc. The position of the 
source can be set as an offset from center, and the angle 
of view can be varied. Models include a Plummer model 
or isothermal sphere. The Einstein radius can be option- 
ally set. A set of default values is also available to show 
the capabilities of the program. 

4.3.3. Summary 

This simulation is simple to install and use. Various 
parameters are easily set in a menu screen and the re- 
sulting lens model instantly visualized. There are two 
optional models that can be used. This software is very 
useful as an educational tool, but not sufficiently robust 
for research analysis. 

4.4. lens 
4.4.1. Software features 

Lens was released in 2002, and is written in MATLAB, 
and therefore the MATLAB (The MathWorks Inc, Nat- 
ick MA USA) software is required. The source code is 
distributed and easily installed in MATLAB (these tests 
conducted with version R2010a). The source code and 
sample data files (SIS1 and SIS2) are available for down- 
load 15 . The package includes two programs written in 
MATLAB, lens ("forward" modeling) and invert (which 
solves the lens inversion problem). The information on 
the website includes a tutorial, and serves as excellent 
documentation. The user interface is through the MAT- 
LAB screen. Installation is accomplished using MAT- 
LAB. Sample data files are provided to test the software. 

15 http : / /www. astro .ubc . ca/people/newbury/siam/lens .html 



11 



1 




\ 


k. 




Fig. 7. — Sample output from lens showing the image plane (left) 
and the generated source plane (right) 

4.4.2. Lensing features 

Lens is a hybrid program, which demonstrates both the 
"forward" problem as well as lens inversion. For model- 
ing the lensed images, lens uses a parametric model as a 
SIS, and the algorithm is described in an accompanying 
article (Newbury 2002) which is easily available follow- 
ing a link from the website. Various model parameters of 
the SIS can be changed in the parameter file, as detailed 
on the website. The input data file specifies observation 
parameters, lens parameters (sigma, z) and source pa- 
rameters (z,x,y, ellipticity, flux) and produces the lensed 
image (the "forward" problem") (left panel, Fig 7). Re- 
sulting images are shown directly on the display with no 
further intervention. The program optionally writes a 
data file which can then be used as input to the inver- 
sion routine to solve the "reverse" problem. The user 
estimates the center of the lensing mass with the com- 
puter mouse, and an image of the source plane is created 
(right panel, Fig 7). 

4.4.3. Summary 

This program is easy to use, but requires MATLAB. 
The program demonstrates a solution to the basic prob- 
lems of gravitational lensing including the forward and 
reverse problems. Input files are easy to create and the 
output graphics are quickly generated and illustrate the 
solution. The software could be easily modified, and is 
supplied with ample comments. 

4.5. MOWGLI 

4.5.1. Software features 

MOWGLI is available as a Java applet, accessed 
through a web browser 16 . The web site includes a Quick- 
start guide on-screen with easy to follow directions. The 
interface is intuitive and requires minimal training to use. 
There is no documentation as all information is described 
on the web site. 

4.5.2. Lensing features 

MOWGLI is designed for interactive manual modeling 
of strong gravitational lenses, as described on the web 
site. Objects such as lenses, sources or masks are added 
to the screen through an easy to use interface. One can 
upload a .jpg image of the system for which modeling 
is desired. By clicking and dragging on the screen, the 
objects can be moved into various positions as desired. 

16 http : //www. ephysics . org/mowgli/ 
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TABLE 2 

Lensing features of gravitational lens modeling software. 

The type of lens models and algorithm used by each of 
the software packages is shown. para=parametric (ltm), 
Non-Para= Non-parametric (Non-LTM) 

Once the objects are positioned, there are two buttons la- 
beled 'predicted' and 'residual' to view these two images 
on the right hand panel of the interface. The residual 
can be minimized by adjusting the Chi square with an 
on-screen button. 

4.5.3. Summary 

This software is easy to use and easily accessible as a 
Java applet. According to the website, a paper has been 
submitted to further describe the features. 

5. DISCUSSION 

Most research using gravitational lens modeling soft- 
ware utilizes a single software package to model an ob- 
served lensing system. Results are then compared with 
observational data or other models from the literature. 
This is the first review to specifically examine strong 
gravitational lens modeling software. This information is 
of particular importance because there are no single stan- 
dard programs used for modeling gravitational lenses. 

5.1. Software Selection 

There are a number of factors which are important 
when selecting appropriate software. Investigators wish- 
ing to develop their own may use existing software, es- 
pecially if source code is provided, as a starting point for 
feature sets, data entry, etc. The computing platform 
may also be an important factor, although most of the 
research software reviewed here will run on a Linux sys- 
tem. Some of the software reviewed runs on a limited 
number of systems. The type of data used for input and 
the results provided are also important factors. Some of 
the programs provide immediate display of output while 
others require another step, with other software, to dis- 
play the results. 

5.2. Software Comparisons 

Comparisons of gravitational lens modeling software 
will be defined as indirect comparisons or direct com- 
parisons. An indirect comparison is comparing re- 
sults from various software, in different papers, mod- 
eling the same lensing system. This is relatively 
easy because the software is commonly tested using 
well-described lensing systems such as C10024+1654 
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(GRALE:(Licsenborgs et al. 2008a)), SDSS J1004+4112 
(glafic:(Oguri 2010), GRALE:(Liesenborgs et al. 2009) 
and PixeLens: (Williams and Saha 2004)), and Abelll689 
(Lenstool: (Julio and Kncib 2009), LensPerfect:(Coe et al. 
2010), and PixeLens: (Saha et al. 2006b)). Direct com- 
parisons compare the models generated for a single lens- 
ing system using different software in the same publica- 
tion. However, there are very few direct comparisons of 
software in the literature. 

The modeling of Abell 1689 by PixeLens, LensPerfect 
and Lenstool has been compared qualitatively (Coe et al. 
2010). Lin and coworkers use Lensview to model a bright 
z=2 galaxy and then compared the results to another 
model made using Lenstool (Lin et al. 2009). There is 
often no single ideal modeling software, as was illustrated 
in the study by Lin and colleagues who found that the 
Lensview model made it difficult to determine how well 
image positions alone are determined (Lin et al. 2009). 
Further studies are needed to compare the ability of var- 
ious software packages to model the same lensing system 
in order to enable meaningful comparisons. Understand- 
ing the underlying assumptions and limitations of the 
various software will be facilitated in the future by more 
direct comparison studies. 

Future research in gravitational lens modeling will also 
require direct comparisons of the results of models using 
data from both strong and weak gravitational lensing 
analysis. Such analyses have already been reported for 
Abelll689 (Limousin et al. 2007; Coe et al. 2010). 

The use of direct comparisons will enable better com- 
parisons of software usability, feature sets as well as 
checking results against known software. Future stud- 
ies in gravitational lens modeling will demand new ap- 
proaches and sophistication. 

6. CONCLUSIONS 

There is a wide variety of gravitational lensing mod- 
eling software available. Many of the publications using 
these packages are written by the software developers, 
suggesting that the software is developed for their per- 
sonal use. Software available only as executable files has 
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the advantage of being rapidly usable, as long as the 
computing platform is available. Software distributed as 
source code may require significant time to compile and 
prepare given the vagaries of software libraries. 

The comprehensiveness of this study is limited by the 
ability to identify existing software or studies related to 
the software reviewed. However, a fairly wide spectrum 
of software was available for review, and this analysis 
identifies opportunities for improvement based on exist- 
ing software. 

Awareness of available software may limit the need to 
develop proprietary software in the future. The source 
code is available for some applications which facilitates 
the development of custom software. Given the increased 
activity in gravitational lensing research, sharing of soft- 
ware and algorithms may result in significant time sav- 
ings. The usability of available software is somewhat lim- 
ited by the essentially consistent use of character based 
user interfaces. Future software should be modular in 
nature and use a graphical user interface for improved 
functionality. 

Until recently, the construction of gravitational mass 
maps to detail dark matter distribution has depended on 
data from weak lensing. However, as shown in recent 
studies using LensPerfect, glafic and GRALE, accurate 
mass maps can be constructed using data from strong 
gravitational lensing data, and this may signal a new 
era in studies of strong gravitational lensing. New ap- 
proaches to software development will be necessary to 
support this shift in research, especially with the advent 
of far more detailed images from the next generation of 
telescopes. Future studies should include direct compar- 
isons with other available software in addition to indirect 
comparisons with previous studies of well-described lens- 
ing systems. 
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